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1 ABSTRACT 

2 Components of the atmospheri c energy budget f rom the M odern Era Retrospecti ve- 

3 anal ysi s for Research and A ppli cations (MERRA) are evaluated in polar regions for the period 

4 1 979-2005 and compared with previous esti mates, i n situ observations, and contemporary 

5 reanalyses. Closure of the energy budget is reflected by the analysis increments term, which 

6 results from virtual enthalpy and latent heating contributions and averages -1 1 W nrf 2 over the 

7 north pol ar cap and -22 W m 2 over the south pol ar cap. T otal energy tendency and energy 

8 convergence terms from MERRA agree closely with previous study for northern high latitudes 

9 but convergence exceeds previ ous esti mates for the south pol ar cap by 46 percent. D i screpanci es 

1 0 with the Southern Hemisphere transport are largest i n autumn and may be related to differences 

11 in topography with earlier reanalyses. For the Arctic, differences between M ERRA and other 

12 sources in TOA and surface radiative fluxes maximize in May. These differences are concurrent 

13 with the largest discrepancies between M ERRA parameterized and observed surface albedo. For 

14 May, i n si tu observati ons of the upwel ling shortwave flux i n the A retie are 80 W m“ 2 larger than 

15 MERRA, whi I e the MERRA downwel ling longwave flux is underestimated by 12 W m“ 2 

1 6 throughout the year. Over grounded ice sheets, the annual mean net surface energy f I ux i n 

17 MERRA iserroneously non-zero. Contemporary reanalyses from the Climate Forecast Center 

1 8 (CFSR) and the I nteri m Re-A nal yses of the European Centre for M edi urn Range Weather 

1 9 Forecasts ( ERA - 1 ) are found to have better surface parameteri zati ons, however these col I ecti ons 

20 are also found to have significant discrepancies with observed surface and TOA energy fluxes. 

21 D i screpanci es among avai I abl e reanal yses underscore the chal I enge of reproduci ng credi bl e 

22 esti mates of the atmospheri c energy budget i n pol ar regi ons. 
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23 1. Introduction 

24 The objective of this study is to examine the performance of theModern Era 

25 Retrospective- analysis for Research and Applications (MERRA) in representing the high latitude 

26 atmospheric energy budget. MERRA was recently released by NASA’s Global Modeling and 

27 Assimilation Office (GMAO). This effort, as well as a companion paper examining the 

28 atmospheric moisture budget (Cullather and Bosilovich, 2010), represent an initial examination 

29 of thi s reanal ysi s i n the pol ar regi ons. 

30 A quantitative knowledge of the flow, storage, and conversion of energy within the 

31 climate system has evolved with time as a result of contributions made by improvements in the 

32 observing system and by numerical atmospheric reanalyses (e.g., Fasullo and T renberth, 2008). 

33 In the polar regions the energy budget and its variability are frequently used as a diagnostic for 

34 understanding rapidly changing conditions including glacial mass balance and perennial sea ice 

35 reducti on (e.g. , Porter et al . , 201 0) . A s noted i n Cul I ather and Bosi I ovi ch (201 0) , numeri cal 

36 reanalyses are widely used in polar research for evaluating polar processes, as boundary 

37 conditions for limited area atmosphere and ocean-sea ice models, and asafirst-order validation 

38 for climate models. However reanalyses inevitably contain inaccuracies resulting from 

39 I i mitations i n the observi ng system, i nconsistencies between differi ng observations, and 

40 incomplete knowledge of the physical processes that are represented in the assimilating weather 

41 forecast model. In particular, surface albedo characteristics over polar oceans and high latitude 

42 cloud properties are both associated with important but complex energy feedback mechanisms 

43 that have historically been poorly simulated (Randall etal., 1998). An initial evaluation of the 

44 high latitude energy budget in areanalysis record is therefore a constructive activity. 

45 Some questions of interest pertaining to this study are as follows. 
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• What are the spati al and temporal patterns of energy budget components i n M ERRA , and 
how do they compare with previous studies and contemporary reanalyses? 

• How do M ERRA surface fluxes compare with in situ field studies? 

• What is the nature of adj ustment terms i n the energy budget? 

Section 2 provides an overview of the M ERRA dataset and method. An evaluation of the 
atmospheric energy balance in polar regions isgiven in section 3. A discussion of these 
comparisons is then given in section 4. 


2. MERRA description and method 

Specification of the MERRA system isgiven in Cul lather and Bosilovich (2010) and 
Rieneckeretal. (2010). MERRA utilizes the incremental analysis update assimilation method 
(IAU; Bloom etal., 1996). In this method, a tendency is computed from the difference between 
an initial 6-hourly analysisfield and the background forecast model state. The forecast model is 
then run a second time over the six-hour interval using this tendency as an additional forcing 
term. The resulting MERRA product isthen composed of dynamical I y-consi stent one-hourly 
fields that are incrementally corrected to observation every six hours. Thus atmospheric budgets- 
as they are constructed in the G EOS-5 AGCM- and their analysis increments are maintained 
within MERRA to the accuracy limited by round-off and data compression errors. 

Following aform similar to Trenberth (1997), the MERRA total energy equation 
integrated over the atmospheric column may be written as 


8 Ae +V-F a =R +F sfc + L v dWv 


dt 


dt 


CHM 


L dw^_ L m 

v dt * dt 


' ANA^e) Qnum ’ 


( 1 ) 


FIL 
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where^/r is total energy in the atmospheric column, is the horizontal transport of total 

68 atmospheric energy, R top isthe downward net radiative flux at the top of the atmosphere (TOA), 

69 F sfc isthe upwelling net surfaceflux, L v isthe latent heat of vaporization, L f is the latent heat of 

70 fusion, W v is column-integrated water vapor (preci pitable water), W t is column-integrated cloud 

71 ice condensate, ANA (E ) isthe sum of analysis increments from the IAU method, and Qnum isthe 

72 sum of spurious numerical adjustments. The time rate of change in total atmospheric energy 

73 storage^ is expressed as 
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( 2 ) 


75 where/tyc is surface pressure, p t0 p is the fixed pressure at the top model level which is 0.01 hPa, 

76 c P is the specific heat of the atmosphere at constant pressure, T v is virtual temperature, <t>s is 

77 surface geopotential , k = \ |v| 2 is ki netic energy, and g is the gravity constant. The product c p T v 

78 is referred to as virtual enthalpy. The divergence term may be expanded as follows: 


79 


Psfc Psfc Psfc Psfc 

V-F,=V- f(Zv<7 v -2y^)V — + V- Lr v V^ + V. fo- V— + V- \k \^~ , (3) 

J g J g J g J g 

Ptop Ptop Ptop Ptop 


80 where <t> is geopotential within the atmospheric column. The net upward surface f I ux is given as 

8 1 =Qh+Qe+ L f p s ~ R sfc ’ ( 4 ) 

82 where Q H and Qe are the upwelling surface turbulent sensible and latent heat fluxes, the product 

83 L f P s is latent heating resulting from solid precipitation, and R s / C isthe net downward radiative 

84 f I ux at the surface. The tendency i mposed by the anal ysi s i ncrements ANA (E) represents the 

85 summation of latent heat, virtual enthalpy, kinetic, and geopotential energy term contributions. 

86 Finally, Qnum'w (1) denotes the contribution of spurious residuals resulting from inertial terms, 
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87 the discretization of the thermodynamic equation, coordinate remapping during model 

88 integration, and time-truncation errors. The relation between M ERRA variables and equation 

89 notation is given in the appendix. 

90 The approach of this study is to evaluate M ERRA values against prior studies for large- 

91 scaleareal averages of thetermsin (1-4) over fixed polar regions as shown in Fig. 1, with a 

92 particular focus on the polar caps. Studies for comparison include Nakamura and Oort (1988), 

93 Genthon and K ri nner ( 1 998) , Serreze et al . (2007) , and Porter et al . (201 0) . N akamura and Oort 

94 (1988) produced budget estimatesfor both polar caps using the ocean flux values of Levitus 

95 (1984), composite satellite data from the period 1966-1977, and atmospheric circulation statistics 

96 from Oort (1983) which are largely based on the upper air station network. Nakamura and Oort 

97 (1988) found the observational network i nsuff icient for computi ng atmospheric energy transport 

98 into the south polar cap and instead produced output from the NOAA Geophysical Fluid 

99 Dynamics Laboratory GCM . Genthon and Krinner (1998) used the 15-year re-analysis of the 

100 European Centre for Medium-Range Weather Forecasts (ECMWF) (ERA-15; Gibson et al., 

101 1997) for the period 1979- 1993 to evaluate the south polar cap. Serreze et al. (2007) examined 

1 02 the north pol ar cap and A rcti c Ocean domai ns usi ng the more recent 40- year re- anal ysi s of 

103 ECMWF (ERA -40; Uppalaet al., 2005) and the National Centers for Environmental 

104 Prediction/National Center for Atmospheric Research Reanalyses (NCEP/NCAR; Kalnay et al., 

105 1996) for the period 1979-2001. Serreze et al. (2007) also examined TOA radiative fluxes from 

106 the Earth Radiation Budget Experiment for the study period February 1985 to April 1989 

107 (ERBE; Barkstrom, 1984). Porter et al. (2010) similarly examined the north polar cap energy 

1 08 budget for the period November 2000 to October 2005 usi ng the 25- year Japanese Re-Analysis 

109 (JRA-25; Onogi et al., 2007), and satellite data from the Clouds and the Earth's Radi ant Energy 
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110 System (CERES; Wielicki etal., 1996) product. In support of budget comparisons with these 

111 previous studies, the evaluation of near-surf ace state variables against in situ station observations 

112 is also instructive. 

113 Corresponding values for surface and TOA energy fluxes are also tabulated for two 

114 contemporary reanalyses for comparison: the ECMWF Interim product (ERA-1; Simmons et al., 

115 2007) and the NCEP Climate Forecast System Reanalysis (CFSR; Sahaet al., 2010). The ERA-1 

116 was produced at T-255 spectral resolution, which is similar to the grid resolution of MERRA. 

117 Energy flux fields are produced from 12- hour forecasts initialized by 4D-Var assimilation. 

118 Monthly fields of the ERA-1 were obtained for the years 1989-2005 at a reduced resolution of 

119 1.5° x 1.5°. The CFSR utilize a coupled atmosphere-ocean model for the initial guessfield and 

120 an interactive sea ice model and was produced at T-382 spectral resolution. Model variables are 

121 produced from 6- hour forecasts. Energy flux fields from the CFSR were obtained at full spatial 

122 resolution. 

123 Evidence of an evolving climate system in polar regions- particularly for the Arctic (e.g., 

124 Porter et al., 2010)- motivate an exclusion of the most recent years in the available MERRA time 

1 25 series for an averagi ng period. The results presented here are for the years 1 979-2005. 

126 

1 27 3. Atmospheric energy budget 

1 28 a. Analysis increments 

1 29 T erms of the atmospheri c energy budget averaged over the peri od 1 979-2005 f rom 

130 MERRA are shown in Table 1 for the polar regions defined in Fig. 1. The far right column 

131 indicates budget adjustment quantities. As noted earlier, artificial moisture filling and chemistry 

132 parameterization terms of the energy budget from equation (1) have essentially zero magnitude. 
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133 Not shown, the spatial pattern of the spurious residual Qmm is characterized by alternating 

1 34 positive and negative val ues i n regions of steep topography. Averages taken over I i mited areas 

1 35 may produce al i asi ng of these osci Mating val ues. For exampl e, Qmm averages 1 .5 W m“ 2 over 

136 the Greenland Ice Sheet. But in general, ANA^ is the largest adjustment quantity of interest in 

137 the atmospheric energy budget, and its spatial patterns are shown in Fig. 2. Here, positive val ues 

138 indicate an energy surplus in the balance equation while negatives indicate a deficit. The 

139 magnitude is a measure of closure obtainable by physical terms. Thespatial patterns shown in 

140 Fig. 2 are complex, vary with time, and are typically dissimilar to the patterns of the analysis 

141 increments for the atmospheric moisture budget shown by Cullather and Bosilovich (2010). 

142 As noted previously, analysis increments for the energy budget are the summation of 

143 contributions from latent heat, virtual enthalpy, kinetic, and geopotential energy terms. Of these 

144 four, the contribution \oANA (E) from virtual enthalpy is largefor monthly and annual averages in 

145 both polar cap regions, while analysis increments from latent heating are also significant for the 

146 north polar cap. For the Northern Hemisphere polar region, negative val ues for ANA (E) arefound 

147 over the A retie Ocean, while positive val ues are present over surrounding lower latitudes. Mean 

148 annual amounts less than -40 W nrT 2 are present in the vicinity of the North Pole with smaller 

149 magnitudes over Greenland and marginal seas. Seasonally, these magnitudes are larger in 

150 summer than in winter, howa/er the val ues do not approach the local imbalances of greater than 

1 51 1 00 W m -2 that are shown for the ERA -40 i n Serreze et al . (2007). For the average over the north 

152 pol ar cap, ANA (E) ranges f rom -4 W m“ 2 i n February to -1 7 W m“ 2 i n June. 

153 In Cullather and Bosilovich (2010), MERRA analysis increments for the atmospheric 

1 54 moisture budget were shown to be characterized by closed contours denoti ng upper ai r stations i n 

155 coastal Greenland and Antarctica Although signatures of upper air station locations are not as 
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156 evident in the ANA (E ) field as in the analysis increments field for the atmospheric moisture 

157 budget, a dipole is apparent in Fig. 2a in the vicinity of Hudson Strait with centers near Kuujjuaq 

158 (58°N, 68°W) and Cape Dorset (64°N, 77°W). Averaged over the north polar cap and the period 

159 1979-2005,^4^; is -11 W m" 2 . The temporal variability of ANA (E ) in the Arctic also differs 

160 markedly from the analysis increments field of the atmospheric moisture budget presented in 

161 Cullather and Bosilovich (2010). Changes in the high latitude atmospheric moisture budget 

162 analysis increments were largely found to be associated with the introduction of datafrom the 

163 Advanced Microwave Sounding Unit (AM SU) in November 1998, which has a global impact on 

164 MERRA (Bosilovich et al., 2010). In contrast, the^A^timeseriesfor the north polar cap 

1 65 i ndicates changes which are not concurrent with satel I ite observi ng system changes. The 

166 magnitude of energy budget analysis increments for the north polar cap averages less than 

167 10W m“ 2 for the period 1 979-1 991 , approxi mately 1 8 W m“ 2 over the period 1 992-1997, and 

168 9 W m“ 2 thereafter. These shifts may be due to changes i n the surface observi ng system or i n 

1 69 atmospheric conditions. 

170 In the Southern Hemisphere, the MERRA energy budget analysis increments term- as 

171 shown in Fig. 2b- hasalarger magnitude than for the north polar cap, with amounts of greater 

172 than (-)80 W rrT 2 over Victoria Land and regions of northern Queen Maud Land in East 

173 Antarctica Over the south polar cap, ANA (E) ranges from -27 W m“ 2 in February to -20 W rrT 2 

174 in August and September. The annual average of ANA (E ) is comparatively smaller over the lower 

175 latitudes of the Southern Ocean asseen in Fig. 2b. There is a considerable annual cycleforthe 

1 76 Southern Ocean domai n rangi ng from -40 W m“ 2 i n January and February to -9 W m“ 2 i n June. 

1 77 The year to year ti me series for the south polar cap is highly variable and ranges from 

1 78 -37 W m -2 i n 1 983 to -9 W m -2 i n 1 998. The anal ysi s i ncrements ti me seri es f or the south pol ar 
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179 cap energy budget is uncorrelated with that of the north polar cap, and its relation to changes in 

1 80 the observi ng system is also not readi ly apparent. But over the data-sparse Southern Ocean 

181 domain there is a discontinuity in \heANA (E ) time series in 1998 that is likely associated with the 

1 82 i ntroducti on of A M SU . Southern Ocean anal ysi s i ncrements average -21 .8 W m“ 2 pri or to 1 998, 

1 83 and -26.7 W rrf 2 thereafter. 

184 

1 85 b. Total atmospheric energy tendency 

186 In both the north and south polar caps, the M ERRA total energy tendency is near zero 

187 for annual averages and is small for months of solstice, as shown in Table 1. But there is an 

188 oscillatory annual cycle for the tendency terms as seen in Fig. 3. For the north polar cap, the 

1 89 tendency term reaches a maxi mum of 26 W nrT 2 i n Apri I and a mi ni mum of -26 W m“ 2 i n 

190 September. This annual cycle agrees very closely with valuesfrom other reanalyses as reported 

191 by Porter et al. (2010), Serrezeet al. (2007), and from the observational study of Nakamura and 

192 Oort (1988). The RMS difference of monthly means with MERRA isonly 4 W m“ 2 for both 

193 NCER/NCAR and JRA-25 as reported by Porter et al. (2010), less than 1 W m -2 for the ERA -40 

1 94 as reported by Serreze et al . (2007), but 1 0 W m“ 2 for Nakamura and Oort (1988). I n general the 

1 95 reanalyses are more si mi I ar to each other than to the earlier Nakamura and Oort ti me series. 

196 For the south polar cap, the total energy tendency in MERRA ranges from a minimum of 

197 -16 W m -2 in April to 30 W m“ 2 in November. Asseen in Fig. 3b, the annual cycle is less 

198 sinusoidal than in the Northern Hemisphere, with the November peak offsetting an average 

199 negative tendency that extends from January through July. The RMS difference with monthly 

200 val ues reported by Nakamura and Oort (1988) as compared to M ERRA is 1 3 W m“ 2 , although 

201 each month is within the standard de/iation of MERRA for the 1979-2005 period. Asseen in 
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202 equation (4), the M ERRA energy tendency incorporates the cloud ice latent heating and kinetic 

203 energy terms which are not considered i n other studies. For monthly means over the regions 

204 examined, however, these terms are negligible. 

205 

206 c. Energy convergence and transport 

207 As seen in equation (3), the divergence term is composed of contributions from latent 

208 heat, virtual enthalpy, kinetic, and geopotential energy terms. For the north polar cap, the annual 

209 cycl e of energy convergence f rom M ERRA consi sts of val ues greater than 1 00 W m“ 2 duri ng 

21 0 wi nter months September through M arch and a mi ni mum of 72 W m“ 2 i n M ay, as seen i n 

21 1 Fig. 3a Porter et al . (201 0) present annual cycles of energy convergence computed as a residual 

212 using several combinations of reanalyses and radiativef lux data sets for the period 2000-2005, 

213 while Serreze et al. (2007) present ERA -40 and NCEP/NCAR reanalysis average monthly values 

214 for the period 1979-2001 . While there is agreement in larger energy convergence in winter, there 

215 is considerable variability among the data sets on the months of the minimum and maximum 

216 value, with May providing aspread of 40 W m“ 2 among the various methods. M ERRA values 

217 concurrent with these previous studies are found within this large range. 

218 Figure 4a indicates that the average poleward energy transport across 70°N iszonally 

219 asymmetric and is focused at preferred longitudes which are associated with the mean longwave 

220 circulation patterns in the middle troposphere (Serreze et al., 2007). In comparison to energy 

221 transports across 70°N from ERA -40 as reported by Serreze et al . (2007) , M ERRA transports 

222 shown in Fig. 4a are comparable but with some differences. First, the poleward (positive) flux 

223 centered near 31 5° E (45° W) has a smal I er zonal extent than i s shown i n Serreze et al . (2007) . 

224 This may be due to the higher spatial resolution of MERRA and the role of the Greenland Ice 
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225 Sheet topography in defining the mid-tropospheric trough pattern over eastern North America 

226 Second, the wintertime poleward transport near 150°E is shown in MERRA to be greater than 

227 20-10 9 W m" 1 . This is stronger by one contour level (5-10 9 W m -1 ) than that shown by 

228 Serrezeet al. for ERA -40. But in general the average meridional transport patterns of MERRA 

229 and ERA -40 are remarkably similar. 

230 I n the Southern Hemisphere, prior studies on atmospheric energy convergence are not as 

231 recent. However comparisons to M ERRA may be made usi ng Nakamura and Oort (1988) and 

232 Genthon and Krinner (1998). Using GCM output, Nakamura and Oort (1988) estimated a mean 

233 annual energy convergence across 70°S of 95 W m“ 2 , which is 23 W m“ 2 less than shown for 

234 MERRA in Table 1. As seen in Fig. 3b, the annual cycle in MERRA contains a broad maximum 

235 over wi nter months and a short period of mi ni mum of val ues less than 1 00 W m“ 2 i n December, 

236 January, and February. In contrast Nakamura and Oort (1988) indicate lower amounts in the 

237 autumn, and their annual cycle is generally more sinusoidal. Nakamura and Oort (1988) and 

238 M ERRA monthly energy convergence values are comparable over the months June to October 

239 but M ERRA is larger by 45 W m“ 2 in January. More recently, Genthon and Krinner (1998) 

240 produced seasonal averages and zonal distributions of energy transport using ERA-15 for the 

241 period 1979-1993. The annual atmospheric energy convergence derived from ERA-15 of 

242 81 W m -2 is considerably smaller than corresponding valuesof either MERRA or Nakamuraand 

243 Oort (1988). Seasonal ly, the largest differences between M ERRA and Genthon and 

244 Krinner (1998) ERA-15 values are in autumn. Energy convergence for the south polar cap for 

245 M arch-Apri I -May averages 134 W m“ 2 in MERRA, while Genthon and Krinner (1998) reported 

246 79W m“ 2 . The spatial distribution of energy transports along the 70°S parallel is strongly 

247 dependent on the meanderi ng coastl i ne, such that spati al resol uti on and topography are 
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248 significant. Thus the differences between M ERRA and ERA-15, though large, may partially 

249 result from differing model grids. Additionally, ERA-15 was known to employ a defective 

250 orography over the ice sheet (Uppalaet al., 2005). A visual inspection of Genthon and Krinner 

251 (1998) results indicates that the ERA-15 mean annual poleward transport is less than M ERRA 

252 near 30° E, an intersection point between the 70°S parallel and the East Antarctic coastal 

253 escarpment. For this location, Genthon and Krinner (1998) plot amounts of between 2 and 

254 3-10 9 W m _1 while M ERRA values are greater than 5-1 0 9 W nf 1 . Additionally Genthon and 

255 Krinner (1998) indicate an annual mean equatorward energy transport in the Ross Sea, while 

256 MERRA indicates an average poleward flux. MERRA and ERA-15 share some general 

257 characteristics of the meridional energy transport including a directional change with season in 

258 the South Pacific region between 180°E and 270°E from poleward during winter months to 

259 equatorward in summer , as shown in Fig. 4b. The figure also shows an opposing seasonal 

260 reversal between 270°E and 300°E in MERRA, and this is also reflected in Genthon and Krinner 

261 (1998). 

262 

263 d. TOA radiative fluxes 

264 For the north polar cap, MERRA TOA radiative fluxes are compared to ERBE (Serreze 

265 et al ., 2007) and CERES (Porter et al ., 201 0). With the exception of midsummer months, the 

266 ArcticTOA radiative flux in MERRA is mainly directed upwards {R top < 0), with an annual 

267 average shown in Table 1 of -110 W m“ 2 . Annual estimatesfrom ERBE and CERES are within 

268 the standard deviation of this value. On the monthly timescale, the largest differences are for the 

269 month of M ay, where the M ERRA 1 979-2005 val ue of -23 W m“ 2 compares with -53 W m“ 2 i n 

270 ERBE (Serreze et al., 2007) and-37W m -2 in CERES (Porter etal., 2010). UsingMERRA 
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271 averages concurrent with these satellite records, MERRA is less than satellite estimates for May 

272 by 29 W m“ 2 as compared to ERBE and by 12 W m -2 as compared to CERES. In July, CERES 

273 indicates a net downwards TOA flux of 21 W m -2 compared to a 1 W m“ 2 upwards flux in 

274 MERRA, while ERBE and MERRA concurrent July values are equal. For other months, the 

275 differences are small. 

276 Table 2 also presents R top values for MERRA in comparison to contemporary reanalyses 

277 of the ERA-1 and CFSR for the period 1989-2005. As seen in Table 2 for the north polar cap, the 

278 MERRA annual net TOA radiative f I ux i s greater than for the other two reanalyses by 4 W m“ 2 . 

279 Again, the largest differences are for the spring time period. For the month of M ay, R top f or the 

280 ERA-I averages —34 W m -2 while the CFSR value is -36 W m -2 , and M ERRA agai n averages 

281 -23 W m“ 2 for 1989-2005. Most of this difference is associated with the upwelling shortwave 

282 flux. For the month of May, the M ERRA TOA upwelling shortwaveflux is less than CFSR and 

283 ERA-1 by 15 W m“ 2 and 18 W m' 2 , respectively. It is noted that CFSR incorporates an 1 1-year 

284 solar cycle while MERRA and ERA-1 use climatological solarforcing, and that ERA-1 usesa 

285 larger solar constant value than M ERRA. For the north polar cap, differences in incoming solar 

286 radiation between M ERRA and CFSR are as large as ±3 W rrT 2 for a given month, but average 

287 lessthanl W m -2 . Differences between ERA-1 andMERRA are as large as 9 W rrT 2 for agiven 

288 month, and average 3 W m“ 2 . 

289 For the south polar cap, theTOA net radiativef lux remains negative throughout the year. 

290 Comparisons to both ERBE data and the values from Nakamura and Oort (1988) indicate that the 

291 annual net TOA radiative flux magnitudein MERRA istoo large, and that the discrepancy is 

292 largest during winter months. The 1979-2005 average net flux as shown in Table 1 for MERRA 

293 is -101 W m -2 . This compares with -90 W m“ 2 from the historical satellite data used in 
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294 Nakamura and Oort (1988), and -95 W m 2 from ERBE for the period February 1985 to April 

295 1989 (Briegleb and Bromwich, 1998). The 1979-2005 annual average for R top is by chance equal 

296 to the 1985- 1989 time peri odf or MERRA. In the annual cycle, the differences with satellite 

297 observations are associated with the winter season. For the months of June, July, and August, the 

298 average flux from Nakamura and Oort (1988) is -131 W m“ 2 , and from ERBE, -134 W m“ 2 . 

299 For M ERRA , the correspondi ng val ue is -1 42 W m“ 2 for both 1 979-2005 and 1 985-1989 ti me 

300 periods. In these winter months, the difference between MERRA and satellite values is almost 

301 enti rely composed of the outgoi ng longwave component. 

302 Values for the south polar cap from ERA -I and CFSR reanalyses tend to agree more 

303 closely with MERRA than with values from satellite data sets. For the 1989-2005 period, the net 

304 TOA radiative flux shown in Table2 averages -101 Wm“ 2 for MERRA, -109 W m -2 for CFSR, 

305 and -102 W m“ 2 for ERA-1. In the annual cycle, the MERRA upwelling longwave flux is less 

306 than the other two reanalyses by 3 W m“ 2 i n spri ng and up to 7 W m“ 2 i n autumn. The CFSR 

307 upwelling shortwave flux is greater than the other two reanalyses by more than 10 W m“ 2 in 

308 December and January, and the ERA-1 solar constant difference with the other two reanalyses 

309 accounts for 3 W m“ 2 on the annual average. 

310 The time series of MERRA TOA radiative fluxes indicate potentially spurious trends in 

31 1 both polar regions. Over north and south polar caps, year to year variability in R top resembles that 

312 of the energy budget analysis increments. For the north polar cap, a maximum for R top is reached 

31 3 ini 993 with -1 07 W m -2 , with val ues as low as -1 1 2 W m“ 2 occurri ng i n 1 981 and 2005. 

314 MERRA TOA fluxesforthesouthpolarcaphaveanirregulartimeserieswitharangebetween 

31 5 mi ni mum and maxi mum val ues of 4 W m“ 2 . Over the Southern Ocean domai n, a sharp change is 

316 noted after 1998. This is likely due to the introduction of AMSU data to the observing system as 
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317 noted earlier. Annual average values prior to 1998 are consistent with an averageof -81 W m“ 2 . 

31 8 For the period 1 999-2005 the M ERRA average for the Southern Ocean is -86 W m“ 2 . 

319 

320 e. Surface fluxes 

321 Figure 5 shows the annual average surface net heat flux from M ERRA for both polar 

322 regions. I n the Northern Hemisphere, smal I negative val ues of between 0 and -5 W irf 2 are 

323 found in a uniform field over nonglaciated land surfaces, which is consistent with subsurface 

324 warming in recent years (Serrezeet al., 2007). Over the central Arctic Ocean, MERRA net 

325 surf ace f I ux val ues are posi ti ve as expected but are excepti onal I y I arge. V al ues greater than 

326 1 5 W m" 2 are found i n the central Arctic, and greater than 20 W m“ 2 i n the approaches to the 

327 North Atlantic. These annual val ues are extraordinary and likely not realistic. A comparison of 

328 the averaged annual timeseriesof monthly val ues with previous studies indicates largest 

329 discrepancies occurring in summer months. The July 1979-2005 net surface flux for MERRA is 

330 -68 W itT 2 for the north polar cap as shown i n Table 1 . This compares with -85 W m“ 2 for 

331 ERA-40 (Serrezeet al., 2007) and -86 W m“ 2 for JRA-25 (Porter etal., 2010). Similar 

332 differences are found between M ERRA and contemporary reanalyses as shown in T able 2. For 

333 the concurrent 1989-2005 averaging period, the July net surface flux for the north polar cap is 

334 -87 W m“ 2 for the CFSR and -78 W m“ 2 for ERA - 1 . 

335 Discrepancies in the surface flux fields are evaluated using observations from the Surface 

336 Heat Budget of the A retie ice camp field study i n the Beaufort Sea i n October 1 997 to October 

337 1998 (SHEBA; Uttal et al., 2002). MERRA surfaceflux val ues are compared with acompilation 

338 of observed SHEBA radiative and turbulent flux measurements by Duynkerke and de Roode 

339 (2001). Comparisons are made using the nearest M ERRA grid point to the reported hourly drift 
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340 camp position. For this data source, SHEBA latent heat flux observations were limited and are 

341 not considered. Usi ng the remai ni ng energy budget components, a net f I ux comparison i ndicates 

342 a positive (upward) biasinMERRA of 18 W m -2 for the months October to April, -1 W m“ 2 for 

343 May, and small positive biasesfor the following summer months. 

344 There are three fundamental results of the comparison. As shown in Fig. 6, substantial 

345 dif f erences i n the upwel ling shortwave radiative flux result from an overly simplistic 

346 representation of sea ice properties. Sea ice albedo is set to afixed value of 0.60 for M ERRA. 

347 The surface observed usi ng SH EBA tower measurements has a much hi gher al bedo i n 

348 springtime, with monthly averages of 0.83 in March, April, and May, and 0.74 in June. Apart 

349 from the tower flux measurements, a line of surface albedo observations made during SHEBA 

350 provide a range of val ues that are dependent on the surface ice conditions. The average of these 

351 surface observations isshown in Fig. 6 for June to September 1998. Albedosfrom tower 

352 measurements i n M ay 1 998 are consistent with area-averaged surface and ai rcraft observations. 

353 For example, Curry et al. (2001) note that albedo for April and May at the SHEBA site averaged 

354 0.84, that the melt season lasted from late May to mid August, and that winter-spring albedo 

355 val ues were again reached in late September. This difference with observed albedo contributes to 

356 an underesti mate i n the upwel I i ng shortwave f I ux i n M ERRA of 55 W m“ 2 i n Apri 1 , 80 W rrT 2 i n 

357 May, and 56 W m“ 2 in June. In late summer, the observed surface albedo isdegraded by melting 

358 and becomes comparable to the M ERRA value. In late autumn, freezing and the introduction of 

359 solid precipitation again produces surface albedo differences between M ERRA and observation, 

360 however the incoming solar flux is reduced and the impact on the upwelling shortwave is less 

361 consequential. The difference with observation in the upwelling shortwave radiative f I ux for 

362 M ay is the largest of any monthly budget component. 
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363 The second result is a response i n other M ERRA surface energy budget terms i n M ay to 

364 the albedo bias. Surface temperatures over icein MERRA are determined via energy balance, 

365 and the underesti mate of the surface albedo results in a perceived increased absorption of solar 

366 energy and a surface warming. This likely results in the springtime MERRA sensible heat flux 

367 bias, whichisfoundtobe16W m -2 in May. Other than the April, May, and June period, the 

368 MERRA sensible heat flux difference with SHEBA observations is only 2 W irf 2 . An intriguing 

369 finding is a springtime negative bias with SHEBA observations in the downwel ling shortwave 

370 radiativef lux. The MERRA downwel ling shortwave is underestimated by 36 W itT 2 in April, 

371 37 W m“ 2 in May and 25 W m -2 in June. In other months this difference is about 1 W m -2 . This 

372 bias is likely associated with general deficiencies in the representation of cloud properties. From 

373 the seasonal timing of the bias, however, it is speculated that a portion of the amount is due to a 

374 redistribution of cloudiness in the atmospheric column resultingfrom anomalous surface 

375 warming. The large May bias in upwelling shortwave radiation is then compensated for by biases 

376 in other fluxes to produce the surface net energy flux bias of -1 W m' 2 . 

377 Shown in Fig. 7 is the time series of hourly near-surface air temperature in comparison to 

378 the observed ti me series from SH EBA for the period 1 February to 30 June, 1 998. A temperature 

379 bias in spring is readily apparent, with a difference of greater than 3.5°C in April and May before 

380 the freezing value is reached in early June. In particular, the period 19 April to 10 May shows an 

381 average bias of 6.1 °C in MERRA. For daily averages, however, there is a good correlation 

382 betweenMERRA and observation for the period shown (r = 0.95). It may be seen in the time 

383 seriesof hourly values shown Fig. 7 that the diurnal cyclein MERRA temperature has an 

384 ampl itude between 2 and 1 0°C, which begi ns abruptly on 28 M arch and conti nues unabated unti I 

385 thefreezing point isreached in June. Theobserved SHEBA diurnal cycle has a similar 
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386 amplitude, however the cycle is not as regular as in MERRA and there are periods of 

387 considerable interruption, perhaps dueto synoptic variability. These differences are suggestive of 

388 difficulties in MERRA boundary layer parameterizations. Springtime air temperature biases are 

389 found at Arctic station locations as well. For example, a comparison with Sachs Harbor (72°N, 

390 125°W) overtheperiod 1979-2005 indicates an average of 4.9°C difference for April but only 

391 1.9°C for the months August through March. A comparison with Barrow (71 °N, 157°W) 

392 similarly indicates an average bias in MERRA of 3.6°C for the spring months of March, April, 

393 and May and 0.9°C for other months. But as shown in Fig. 8, MERRA performs well in a 

394 comparison of monthly anomalies. The correlation between temperature anomalies at Barrow 

395 and Jan Mayan (71°N, 9°W) is0.99for both stations. The time series shown in Fig. 8 contain 

396 observations that cover the entire time period. Other stations in the Arctic with shorter and/or 

397 i nterrupted records compare similarly well. 

398 The third result from the comparison with SHEBA is a negative bias in the downwel ling 

399 longwave radi ative f I ux throughout the year of 1 2 W m' 2 . This quantity leads to the overal I 

400 positive bias in the net surface flux for summer, autumn, and winter months. As with the 

401 springtime downwel ling shortwave radi ative flux bias, an inadequate representation of cloud 

402 properties is implied. To evaluate this further, comparisons were made between MERRA and 

403 SHEBA hourly microwave radiometer retrievals over the period 5 December 1997 to 

404 9 September 1 998. M ore than 5000 observati ons were made over the peri od. Retri eval s of 

405 precipitable water compare remarkably well to MERRA values as seen in Fig. 9a, although 

406 differences are apparent for small quantities in winter. For monthly intervals, the correlation 

407 between MERRA and the hourly microwave radiometer precipitable water retrievals ranges from 

408 r = 0.87 in December 1997 to r=0.96 in May 1998. A consistent bias of 0.6 mm in monthly 
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409 averages is found, which amounts to 31 percent of the observed average for January but only 

410 3 percent for July. In contrast, the comparison to retrieved liquid water content shown in Fig. 9b 

411 is less favorable. Cloud liquid water from the SHEBA microwave radiometer rangesfrom an 

412 average of 0.017 mm in January 1997 to 0.106 mm in August 1998. Typical MERFtA values are 

413 about 45 percent of the microwave radiometer amounts. Although large discrepancies have been 

414 noted between the SHEBA microwave radiometer valuesfor liquid water path and simultaneous 

415 aircraft measurements (Lin et al., 2001), the differences between MERFtA and SHEBA values 

416 exceed 50 percent. Additionally, the correlations of hourly liquid water path values with 

417 MERRA over monthly time intervals are low and range from r = 0.1 4 in April 1998 to r= 0.55 in 

418 January 1998. The presence or absence of cloud liquid water significantly alters the downwel ling 

419 longwave radiative flux. An underestimate of cloud liquid water in MERRA is qualitatively 

420 consistent with differences in the surface net flux with observation. 

421 Comparisons with MERRA for the Arctic are also conducted using the CFSR and ERA-1 

422 reanalyses. Using monthly values co-located with the SHEBA ice drift camp, it is noted that 

423 surface albedo varies seasonally and interannually in both CFSR and ERA-1. In agreement with 

424 the SHEBA time series, both CFSR and ERA -I have albedos greater than 0.8 for April 1998, and 

425 values decrease with the onset of the summer melt season. This decrease occurs more rapidly in 

426 both CFSR and ERA-1 than for tower observations, but is within the lower range given by 

427 SHEBA line albedo measurements. The June 1998 albedo is 0.59 for MERRA, 0.65 for CFSR, 

428 0.69 for ERA -I, 0.74 for the SHEBA tower observation, and 0.62 for the line observation. All 

429 three reanalyses underestimate the downwel ling longwave radiative flux overwinter months in 

430 comparison to SH EBA . For the period October 1 997 to M ay 1 998, this f I ux is underesti mated by 

431 5 W m 2 in ERA-1 and 18 W m 2 in CFSR. Finally, the November 1997 to March 1998 average 
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432 sensible heat flux observed at SHEBA is less than 1 W m“ 2 . This compares with 3 W m“ 2 in 

433 MERRA, -7 Wm" 2 in ERA -I, and -2 1 W m“ 2 i n CFSR. 

434 T urni ng to the Southern Hemisphere, the annual average net surface heat f I ux for the 

435 south polar cap is shown in Fig. 5b. Of immediate concern is the anomalous non-zero field over 

436 Antarctica, which is shared by the major ice sheets in both polar regions. Over grounded ice, the 

437 M ERRA subsurface energy f I ux is determi ned i n the GEOS-5 model by the prognostic 

438 temperature for a 7 cm (water-equivalent) surface ice layer and a “deep” layer temperature at 2 m 

439 depth that isfixed at 230°K. Thus, the location of the zero value contour in Fig. 5b exactly 

440 matches the annual -average 230°K surface temperature isotherm. Observations from automatic 

441 weather stations indicate that annual mean subsurface conductive heat fluxes are not significant 

442 (e.g., Reijmer and Oerlemans, 2002), and annual surface energy flux patterns in M ERRA over 

443 Antactica (as well as Greenland) are erroneous. 

444 The pattern in the MERRA annual surface net energy flux in Antarctica is manifest as a 

445 complementary distribution of downward (negative) turbulent and upward (positive) radiative 

446 fluxes that are not balanced. MERRA annual mean latent heating exceeds 5 W m“ 2 only along 

447 the East Antarctic coast i n selected locations, and averages less than 1 W m“ 2 for the total 

448 grounded ice sheet area The annual averaged sensible heat flux over the ice sheet is uniformly 

449 negative and is approximately contour-parallel with topography, with magnitudes greater than 

450 (-)60 W itT 2 al ong the East A ntarcti c coastal escarpment decreasi ng to I ess than (-) 1 0 W rrT 2 

451 over the central plateau. The annual mean net radiative flux field in MERRA is spatially more 

452 uniform with values ranging from (+)25 to 35 W m“ 2 for East Antarctica, and smaller positive 

453 val ues over West A ntarctica This results i n the i mbalances i n the net surface heat f I ux as shown. 

454 Over the interior plateau, net f I ux val ues are as large as +15 W m“ 2 while the net amounts at 
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455 lower elevations are negative and are less than -30 W m“ 2 over the East Antarctic coastal 

456 escarpment. These errors in the net surface flux have relation to near-surface temperature biases. 

457 Asshown in Fig. 10, there is a considerable wintertime warm bias of 5°C at Amundsen-Scott 

458 (90°S), while a summer cold bias of 5°C is found at Scott Base(78°S, 167°E). Visual 

459 comparison with satellite-derived surface air temperatures in Comiso (2000) indicates that a 

460 summer cold bias extends over the embay ment regions. 

461 Comparisons of surface energy budget components are made with Antarctic station 

462 values compiled by King and Turner (1997). Values compiled by King and Turner (1997) reflect 

463 studies of opportunity and do not account for interannual variability. For the sensible heat flux, 

464 M ERRA averages at M izuho (71 °S, 44°E) of -47 W m“ 2 i n July and -1 9 W m“ 2 i n December 

465 compare with observational values of -37 and -25 W rrf 2 for July and December, respectively 

466 (Ohataet al., 1985). At South Pole, differences in seasonal values of the sensible heat flux are 

467 largest in spring and summer. The December-January- February sensibleflux averagefrom 

468 MERRA is-3W m“ 2 and-22W m“ 2 in observation (Carroll, 1982). This contributes to a 

469 difference of 9 W m -2 in the annual average. 

470 Differences between M ERRA and available observations are also associated with the net 

471 radiative flux. At Halley (76°S, 26°W), the annual average net radiativef lux for MERRA of 

472 13 W rrf 2 approximates the observational values of 9.8 W m“ 2 , however seasonal differences are 

473 as large as 10 W m -2 in winter. At South Pole, the annual net radiative flux of 19 W m“ 2 matches 

474 the observation of Carroll (1982), however seasonal differences are large. In winter, the net 

475 radi ati ve surface cool i ng of 36 W m“ 2 exceeds the observed val ue of 21 W m“ 2 . 1 n summer, the 

476 MERRA radi ative flux value of -7 W m“ 2 differs from the Carrol I (1982) valueof +18W m“ 2 . 
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477 The surface radiativeflux differences at South Pole are further examined using the 

478 observati ons of D utton et al . ( 1 989) , w ho recorded dai I y mean radiative flux components f rom 

479 April 1986 until February 1988. Over this period, the surf ace net radiativeflux for both MERRA 

480 and observation is positive for most of the year but becomes negative i n summer months as seen 

481 in Fig. 11. Over the 22 month period, the downwel ling longwave radiativeflux is consistently 

482 less than observation by an average 24 W m“ 2 . This difference is apparent i n the comparison of 

483 daily valuesin Fig. 11. Large biases are also found intheMERRA net shortwave radiativeflux 

484 i n spri ng and summer. For the month of January, the net downward shortwave f I ux is 

485 overesti mated by 20 W m 2 i n 1 987, and by 23 W m“ 2 i n 1 988. A mi nor part of the shortwave 

486 bias is associated with the MERRA surface albedo, which isfixedover land ice at 0.775. 

487 Observed monthl y averages at South Pol e i ndi cate an al bedo of between 0.80 and 0.89. These 

488 differences in the shortwave flux partially cancel the downwel ling longwave underestimate in 

489 summer. It may be seen from Fig. 11 thatsomeof the day-to-day variability in the downwel ling 

490 longwave radi ative f I ux is reproduced i n M ERRA . By subtracti ng a 30-day runni ng mean from 

491 each ti me seri es to remove the annual cycl e, the two curves have a correl ati on of 0.70. 

492 Table 2 presents a comparison of net surf ace flux values for the south polar cap. Boththe 

493 ERA-1 and the CFSR correctly depict a near-zero annual net flux field over the Antarctic ice 

494 sheet, while regions of opposite sign in M ERRA Fsfc fortuitously cancel. Monthly values of 

495 surface radi ative flux components from CFSR and ERA-1 are compared to 1986- 1988 values 

496 from Dutton et al. (1989) for the South Pole. The ERA-1 collection begins in 1989, so 1989-2005 

497 averages were used. In general, the monthly net surface radi ative fluxes of the three reanalyses 

498 are more similar to each other than to observation. The net upward radiativeflux is 

499 overestimated by 18 W m -2 for MERRA, 16 W m“ 2 for ERA-1, and 20 W m“ 2 for CFSR. Smilar 
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500 to MERRA, alargepart of the ERA-1 difference is due to an underestimate of the downwel ling 

501 longwave component. For the annual average, the ERA -I downwel ling longwave flux is 

502 underestimated by 15 W m“ 2 . FortheCFSR, theupwelling longwaveflux is overestimated for 

503 winter months March to September by 21 W m -2 , and this provides a significant contribution to 

504 annual net flux differences. Both ERA-1 and CFSR have variable surface albedos at South Pole, 

505 however the ERA -I value approximates the MERRA fixed value, while the CFSR value averages 

506 0.84 duri ng summer months. 

507 The spati al patterns of Fi g. 5b are of i nterest over the Southern Ocean. I n the annual 

508 mean, M ERRA i ndicates a net loss of energy from the ocean to the atmosphere south of 60°S 

509 which increases in magnitude near the continent. Farther north there is a marked asymmetry 

51 0 withi n the 50°S-60°S zone, with net energy loss from the ocean to the atmosphere i n the Pacific 

511 sector and energy gains elsewhere. Embedded within the Pacific sector are two regions of net 

512 energy gain from the atmosphere to the ocean which correspond to meanderings of the Antarctic 

513 Polar Front- as it crosses the Southeast Indian Ridge near 145°E, and the Pacific- Antarctic Ridge 

514 near 145°W (e.g., Moore et al., 1999). Josey (2009) noted the zonal asymmetry in the net surface 

515 heat flux in NCEPand ECMWF reanalyses but found that coupled models produce a more 

516 zonally uniform field. Josey (2009) concluded that the sign of annual mean surface heat 

517 exchangeover much of the region is not known. The net surfaceflux from ERA-1 and CFSR 

518 reanalysesin Southern high latitudes differ with MERRA. Thecoastal zoneof heat lossfromthe 

519 ocean to the atmosphere i n both the ERA-1 and CFSR is more closely confined near the continent 

520 than in MERRA. Similar to MERRA, the ERA -I indicates an annual mean net positive energy 

521 f I ux from the ocean to the atmosphere i n the Pacific Ocean sector of the 50°S-60°S zone, whi le 
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522 theCFSR indicates negative values between 0 and -15 W nT 2 that are smaller in magnitude than 

523 for the rest of the zone. 

524 The annual cycle of the net surface f I ux for the Southern Ocean is shown i n Fig. 3c. 

525 Using ECM WF operational analyses overlapping the period of the ERBE study, Okada and 

526 Yamanouchi (2002) examined the atmospheric energy budget for the region bounded by 60°S 

527 and 70°S. Okada and Y amanouchi esti mated the surface energy budget as the residual usi ng 

528 TOA ERBE radiation and analyses divergence terms. A seasonal asymmetry in the net surface 

529 flux was highlighted, which was found to abruptly peak in May with a maxi mum value of 

530 1 1 6 W rrf 2 . Okada and Y amanouchi (2002) attributed this asymmetry to the latent heat release 

531 resulting from sea ice formation. As seen in Fig. 3c, theMERRA surface energy flux over the 

532 Southern Ocean sea i ce domai n i s al so asymmetri c and peaks i n M ay at 98 W m“ 2 , however the 

533 maximum is not as striking as was found for the ECM WF analyses. In examining the autumnal 

534 surface turbulent fluxes in MERRA, it is found that the total latent heatflux isamaximumfor 

535 the domain in April with 33 W m“ 2 . The latent heat flux then diminished over ice covered winter 

536 months, with a second maxi main November of 28 W m“ 2 . TheMERRA sensible heat flux 

537 reaches its annual maximum in May of 21 W m“ 2 and generally reflects the shape of surface net 

538 flux. The asymmetry in the annual cycle for the MERRA net surfaceflux asshown in Fig. 3c is 

539 principally due to seasonal changes in the sensible heat flux. In reanalyses, sea ice cover is 

540 prescribed from observational fields. The latent heatflux arising from iceformation is manifest 

541 as the net conductive flux at the atmosphere- ice interface. In this context, MERRA and the 

542 results of Okada and Yamanouchi (2002) are broadly consistent. 

543 
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544 4. Summary and Discussion 

545 M ERRA reproduces the basi c patterns of energy f I ow i n the pol ar atmosphere as they are 

546 known. Asshown in Fig. 3, the polar regions are marked by a convergence of energy from lower 

547 latitudes for all months, and a loss of energy at the top of the atmosphere for the most of the year. 

548 IntheArctic, reductions in the TOA shortwave radiative flux in autumn produce a negative 

549 tendency in the atmospheric column total energy throughout the period August through January, 

550 which is moderated by contributions from the net surface flux and increased energy transport 

551 from lower latitudes in winter (Serrezeet al., 2007). In the Antarctic, this seasonal progression is 

552 I ess sinusoidal, with the net TOA radiative flux remaining negative throughout the year, and an 

553 extended wi nter peri od i n the energy budget components extendi ng f rom A pri I to September. 

554 Despite reproducing these essential components, MERRA energy budgets for the Arctic 

555 and A ntarcti c contai n substanti al errors owi ng to overl y simplistic physi cal parameteri zati ons, 

556 including sea ice albedo, the surface heat budget over permanent land ice, and cloud radiative 

557 properties. Difficulties in MERRA with sea ice characteristics are not dissimilar from those 

558 descri bed i n Bretherton et al . (2000) for ECM WF anal yses produced duri ng SH EBA , and i ndeed 

559 the discrepancies in surface shortwave radiative fluxes are similar. Spring is acritical period for 

560 evaluation of surfaceflux fields in the Arctic, and differences between MERRA shortwave 

561 surface radi ative f I uxes with observation are most promi nent i n M ay. Over the data sparse 

562 Southern Ocean, discontinuities in the time series of TOA radi ative fluxes coincide with the 

563 introduction of AMSU satellite data in November 1998 and are therefore spurious. Elsewhere, 

564 interannual variability of the analysis increments term ANA (E ) is large but not as easily linked to 

565 changes in the observing system. Additional characterization of analysis increments, including 
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566 their vertical distribution, and MERRA cloud properties are conspicuous points for further 

567 evaluation. 

568 M ERRA nevertheless compares favorably to previous studies of energy budget 

569 components produced from state and dynamical variables. These vertical integrals are pre- 

570 computed quantities in MERRA, and are not readily aval I able from contemporary reanalyses. 

571 Atmospheric energy convergence and the spatial distribution of transport along the 70° parallel 

572 compare closely with previous studies in the Northern Hemisphere, while estimates for the south 

573 polar cap are qualitatively similar but may also been seen as an update to studies based on earlier 

574 analyses. The total atmospheric energy tendency in polar regions also compares favorably to 

575 previous studies. 

576 Credible estimates of the atmospheric energy budget in polar regions continue to bea 

577 significant challenge due to changes in the observing system and complex energy feedback 

578 mechanisms that are associated with the high latitudes. Evaluation using both representative 

579 point location observations and previous area-averaged estimates such as those used in this study 

580 are valuable for providing a straightforward appraisal of new reanalyses. The M ERRA system is 

581 an important product due to its alternative construction, including a non-spectral background 

582 model and its emphasis on NASA satellite products. An important concept used in M ERRA is 

583 the employment of analysis incrementsfor identifying differences between observations and the 

584 background analysis system. Inconsistencies in atmospheric budgets are quantified in the 

585 analysis increments, which is one means of measuring uncertainty. ERA-1 and CFSR reanalyses 

586 arefoundto utilize seasonal variations in sea ice albedo and have realistic annual mean surface 

587 heat fluxes over ice sheets. However these col lections are also found to have significant 

588 discrepancies with observed surface and TOA energy fluxes. In particular, sensible heat fluxes 
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589 from CFSR are large in comparison to SHEBA observations, while all three reanalyses 

590 overestimate the annual surface net radiative flux at South Pole by 16 to 20 W m“ 2 . These 

591 disagreements underscore the challenge of the high latitude energy budget problem. 

592 A general criticism of reanalyses is that they are produced with the i ntent of providi ng the 

593 best representation of conditions for a given time without consideration for the impact of 

594 heterogeneous observations on temporal variability (Thorne and Vose, 201 0). This intent 

595 nevertheless has practical, scientific application. Additionally, reanalyses may beseen as part of 

596 a spectrum of products for climate study ranging from heterogeneous observations to model 

597 simulations, which include AM IP fields and sparse data reanalyses (e.g., Compoet al., 2006). As 

598 part of that continuum, the analysis increments in M ERRA provides a quantification of 

599 differences between observations and the background system. Changes in the spatial and 

600 temporal variability of the analysis increments imply changes to the observing system, which 

601 should be carefully treated in evaluating time series. MERRA is avaluable record for examining 

602 the polar atmosphere when these cautions are exercised. 

603 
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616 APPENDIX 

617 


618 

619 

620 
621 
622 

Representation of the Atmospheric Energy Budget Using MERRA Variables 
Thefollowing MERRA variables are given asfollows. 

DQVDT DYN Vertically integrated water vapor tendency for dynamics 

kg rrf 2 s' 1 

623 

DQVDTPHY 

Vertically integrated water vapor tendency for physics 

kg rrf 2 s' 1 

624 

DQVDTANA 

Vertically integrated water vapor tendency for analysis 

kg rrf 2 s' 1 

625 

DQIDTDYN 

Vertically integrated ice water tendency for dynamics 

kg rrf 2 s' 1 

626 

DQIDTPHY 

Vertically integrated ice water tendency for physics 

kg rrf 2 s' 1 

627 

DQIDTANA 

Vertically integrated ice water tendency for analysis 

kg rrf 2 s _1 

628 

DQVDTCHM 

Vertically integrated water tendency for chemistry 

kg m“ 2 s _1 

629 

DQVDTFIL 

Artificial “filling” of water vapor 

kg m“ 2 s _1 

630 

DQIDTFIL 

Artificial “filling” of frozen water 

kg m“ 2 s _1 

631 

DKDTDYN 

Vertically integrated kinetic energy tendency for dynamics 

WrrT 2 

632 

DKDTPHY 

Vertically integrated kinetic energy tendency for physics 

WrrT 2 

633 

DKDTANA 

Vertically integrated kinetic energy tendency for analysis 

WrrT 2 

634 

DHDTDYN 

Vertically integrated c P -T v tendency for dynamics 

WrrT 2 

635 

DHDTPHY 

Vertically integrated c P -T v tendency for physics 

WrrT 2 

636 

DHDTANA 

Vertically integrated c P -T v tendency for analysis 

WrrT 2 

637 

DPDTDYN 

Potential energy tendency for dynamics 

WrrT 2 

638 

DPDTPHY 

Potential energy tendency for physics 

WrrT 2 

639 

DPDTANA 

Potential energy tendency for analysis 

WrrT 2 

640 

CONVKE 

Vertically integrated convergence of kinetic energy 

WrrT 2 

641 

CONVCPT 

Vertically integrated convergence of virtual enthalpy 

WrrT 2 

642 

CONVPHI 

Vertically integrated convergence of geopotential 

WrrT 2 

643 

SWTNT 

TOA outgoing shortwave flux 

WrrT 2 

644 

SWGNT 

Surface net downward shortwave flux 

WrrT 2 

645 

LWTUP 

Upward TOA longwaveflux 

WrrT 2 

646 

LWGNT 

Net downward longwaveflux at the surface 

WrrT 2 

647 

EFLUX 

Latent heat flux (positive upward) 

WrrT 2 

648 

HFLUX 

Sensible heat flux (positive upward) 

WrrT 2 

649 

PRECSN 

Frozen precipitation at the surface 

kg m“ 2 s _1 
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Generation of kinetic energy 

Total energy added by artificial energy “fixer” 


650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680 
681 
682 


DKDTGEN 

TEFIXER 


W m“ 2 
WitT 2 


A tendency may be expressed as the sum of dynamics, physics, and analysis variables. For 
example, the tendency of vertically integrated water vapor (preci pi table water) is expressed using 
MERRA variables as follows. 

dW 

— := DQVDTDYN+DQVDTPHY+DQVDTANA (5) 

dt 


Equation (2) is represented as follows. 

c) 

— ^ := L v • ( DQVDTDYN+ DQVDTPHY + DQVDTANA) 

dt 


i L r {DQIDTDYN + DQIDT PHY + DQIDTANA) 
+ DHDTDYN + DHDTPHY + DHDTANA 
+ DPDTDYN+ DPDTPHY+ DPDTANA 
+ DKDT DYN+ DKDT PHY + DKDT ANA 


Equation (3) is represented as: 

V • F„ :=-(L v • DQVDT DYN - L r DQIDT DYN 
+ CONVKE + CONVCPT+ CONVPHI) 

Equations (4) and (5) are represented as follows. 

Rtop + F sfc := {SWTNT - LWTUP) i (SWGNT + LWGNT) 
+ EFLUX+ HFLUX+ L r PRECSN 

The contribution of spurious residuals in the energy term is represented as: 

Qnum := i DKDT DYN + CONVKE + CONVPHI 
+ DKDT GEN i DPDT DYN i TEFIXER 


( 6 ) 


(7) 


( 8 ) 


(9) 


The remainder of equation (1) is given as follows. 


8W„ 


dt 


+ 

CHM L 


L^-L 
v dt h 


dW t 

dt 


+ ANA {e) \- 


J l FIL 


L v • DQVDT CHM+ (L v • DQVDT FIL - L r DQIDT FIL) 
+ (L v • DQVDT_ANA i L r DQIDT ANA + DHDT ANA 
+ DKDTANA + DPDTANA ) 


( 10 ) 
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800 Table 1 . Components of the M ERRA atmospheric energy budget for regions defined in Fig. 1 , 

801 in W m“ 2 . The surface flux F sfc discounts latent heating from solid precipitation. The standard 

802 devi ati on over the 1 979-2005 ti me peri od i s i ndi cated i n parentheses. 

803 

804 Table2. MERRA, CFSR, and ERA-1 1989-2005 averageTOA and surface energy flux values 

805 for regions defined in Fig. 1, in W m“ 2 . The standard devi ati on over the time period is indicated 

806 in parentheses. 

807 

808 List of Figure Captions 

809 Figure 1 . Regions of study for (a) the Northern Hemisphere and (b.) the Southern Hemisphere. 

810 Bold line indicates the 70° parallel. Continental areas are shaded gray. 

811 

812 Figure 2. Average MERRA analysis increments field for the atmospheric energy budget 

813 (variable>4jV^) for (a) the Northern Hemisphere and (b.) the Southern Hemisphere. The 

814 contour interval is 20 W m -2 . The zero contour is indicated with asolid black line. 

815 

816 Figure 3. Annual cycleof atmospheric energy budget components in MERRA for (a) north 

817 polar cap, (b.) south polar cap, and (c.) the Southern Ocean domain, in W m“ 2 . Bars indicate plus 

818 and minus the standard deviation for the period 1979-2005. 

819 
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820 Figure 4. Average monthly meridional energy transport from M ERRA (a) across 70°N, 

821 contoured every 5-1 0 9 W rrf 1 , and (b.) 70°S, contoured every 3*1 0 9 W nrf 1 . Positive values 

822 i ndicate northward transport. 

823 

824 Figure 5. Annual average net surface heat flux from M ERRA (positive upwards). Contours are 

825 plotted with an interval of 20 W m" 2 and for the levels -10, -5, 0, 5, and 10 W m“ 2 . The zero 

826 contour is indicated with a sol id black line. 

827 

828 Figure 6. Monthly averaged surface albedo (gray) and upwel ling shortwave radiative flux (dark) 

829 for SH EBA observed (sol id) and correspondi ng M ERRA val ues (dashed) for October 1 997 to 

830 September 1 998, i n W m" 2 . “Tower” values are from downward-poi nti ng pyranometer 

831 measurements, while “line albedo” val ues are from surface measurements along a 300m I i ne. 

832 

833 Figure 7. Near-surface hourly air temperature from SHEBA and corresponding valuesfrom 

834 M ERRA for the period 1 February 1 998 to 30 June 1 998, i n degrees C. 

835 

836 Figure 8. Time series of monthly averaged near-surface station air temperature anomaly and 

837 corresponding M ERRA valuesfor Barrow (left, 71°N, 157°W) and Jan Mayen (right, 71°N, 

838 9° W), in degrees C. 

839 

840 Figure 9. Hourly (a) preci pi table water and (b.) liquid water path from SHEBA microwave 

841 radiometer and corresponding M ERRA values, in mm. 

842 
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843 Figure 10. Average annual time series for near surface station temperature and corresponding 

844 M ERRA values for (a) Amundsen-Scott (90°S), and (b.) Scott Base (78°S, 167°E), in degrees C. 

845 Bars indicate the standard deviation of monthly values over the period 1979-2005. 

846 

847 Figure 11. Time series of daily downwelling longwave flux and the net downward flux from 

848 Dutton et al . (1989) and correspondi ng val ues from M ERRA for 90°S, i n W m“ 2 . 

849 
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Table 1 . Components of the M ERRA atmospheric energy budget for regions defined in Fig. 1 , 
in W m -2 . The surface flux F sfc discounts latent heating from solid precipitation. The standard 
deviation over the 1979-2005 time period is indicated in parentheses. 





-V'F a 

Rtop 

Fsfc 

ANA(e) i Qmjm 

70°N -90°N 

January 

-1(11) 

110(16) 

— 1 73( 4) 

63( 5) 

-5( 9) 


July 

2( 6) 

81 ( 8) 

2( 3) 

— 68( 4) 

-13( 5) 


Mean 

0( 2) 

99( 4) 

— 1 1 0( 1) 

19( 1) 

-11( 5) 

Arctic Ocean 

January 

-3(13) 

106(19) 

— 1 76( 4) 

70( 5) 

-5(14) 


July 

3( 9) 

96(10) 

-7( 4) 

-75( 6) 

-9( 6) 


Mean 

0( 2) 

99( 6) 

— H4( 1) 

23( 2) 

-9( 7) 

Greenland 

January 

-2(19) 

143(26) 

— 1 53( 6) 

-1( 4) 

4(16) 


July 

1( 8) 

118(19) 

— 48( 1) 

-40 ( 1) 

-32(16) 


Mean 

0( 2) 

138( 9) 

— 1 1 2( 1) 

— 1 6( 1) 

1 4( 9) 

70°S-90°S 

January 

-1(10) 

78(10) 

-25( 2) 

-32( 2) 

-23( 8) 


July 

-10(13) 

131(16) 

-142( 3) 

19( 2) 

-21(10) 


Mean 

0( i) 

1 1 8( 6) 

— 1 01 ( 1) 

3( 1) 

-22( 6) 

Southern 

January 

2( 6) 

97(17) 

58(10) 

-114( 9) 

-40(17) 

Ocean 

July 

-7( 8) 

91(15) 

-174( 1) 

82( 3) 

-11(13) 


Mean 

0( 1) 

89( 5) 

-83( 2) 

12( 2) 

— 23( 5) 

Antarctica 

January 

-2(10) 

80(11) 

-41 ( 1) 

-21 ( 1) 

-21 ( 8) 


July 

-9(12) 

135(18) 

— 134( 3) 

7( 2) 

-20(10) 


Mean 

0( 1) 

124( 5) 

— 1 01 ( 1) 

-3( D 

— 23( 5) 




Table2. MERRA, CFSR, and ERA-1 1989-2005 averageTOA and surface energy flux values 
for regions defined in Fig. 1, in W m~ 2 . The standard deviation over the time period is indicated 
in parentheses. 





Rtop 



F s fc 


MERRA 

CFSR 

ERA-1 

MERRA 

CFSR 

ERA-1 

70°N -90°N 

January 

-172(4) 

-172(5) 

-174(4) 

64(6) 

55(7) 

59(7) 


July 

2(3) 

-1(3) 

0(3) 

-69(4) 

-87(5) 

-78(4) 


Mean 

-110(2) 

-114(1) 

-114(1) 

19(1) 

14(1) 

12(2) 

Arctic Ocean 

January 

-176(4) 

-174(5) 

-177(4) 

70(5) 

59(7) 

55(5) 


July 

-5(4) 

-2(3) 

-1(3) 

-79(6) 

-106(8) 

-90(6) 


Mean 

-114(1) 

-115(1) 

-117(1) 

23(1) 

12(2) 

12(2) 

Greenland 

January 

-152(5) 

-159(5) 

-159(5) 

0(3) 

6(1) 

11(2) 


July 

-48(1) 

-58(2) 

-51(2) 

-40(1) 

-12(2) 

-15(1) 


Mean 

-112(1) 

-120(2) 

-116(2) 

-16(1) 

1(0) 

2(0) 

70°S-90°S 

January 

-26(2) 

-34(3) 

-20(2) 

-31(2) 

-19(2) 

-22(2) 


July 

-142(2) 

-147(3) 

-146(3) 

19(2) 

16(4) 

12(2) 


Mean 

-101(1) 

-109(1) 

-102(1) 

3(1) 

7(2) 

3(1) 

Southern 

January 

55(11) 

76(5) 

72(3) 

-111(9) 

-145(5) 

-133(3) 

Ocean 

July 

-173(1) 

-172(1) 

-177(1) 

81(10) 

84(5) 

67(3) 


Mean 

-84(3) 

-77(2) 

-84(1) 

13(2) 

3(5) 

-1(1) 

Antarctica 

January 

-41(1) 

-54(2) 

-38(2) 

-21(1) 

-2(0) 

-5(1) 


July 

-134(3) 

-140(3) 

-138(3) 

7(2) 

2(0) 

3(1) 


Mean 

-101(1) 

-109(1) 

-100(1) 

-3(1) 

1(0) 

0(0) 
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Figure 2. Average M ERRA analysis increments field for the atmospheric energy budget 
(variable>4iV^ r£ ;) for (a) the Northern Hemisphere and (b.) the Southern Hemisphere. The 
contour interval is 20 W m -2 . The zero contour is indicated with asolid black line. 
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Figure 3. Annual cycle of atmospheric energy budget components in MERRA for (a) north 
polar cap, (b.) south polar cap, and (c.) the Southern Ocean domain, in W rrT 2 . Bars indicate plus 
and minus the standard deviation for the period 1979-2005. 



41 





Month 



. 


24 

21 

18 

15 

12 

9 

6 

3 

0 

-3 

-6 

-9 

-12 

-15 

-18 

-21 

-24 


70°S (b.) 

Figure 4. Average monthly meridional energy transport from MERRA (a) across 70°N, 
contoured every 5*1 0 9 W m 1 , and (b.) 70°S, contoured every 3*1 0 9 W m 1 . Positive values 
indicate northward transport. 


42 



100 

80 

60 

40 

20 

15 

10 

5 

0 

-5 

-10 

-15 

-20 

-40 

-60 

-80 

-100 



Figure 5. Annual average net surface heat flux from M ERRA (positive upwards). Contours are 
plotted with an interval of 20 W m“ 2 and for the levels -10, -5, 0, 5, and 10 W m“ 2 . The zero 
contour is indicated with a sol id black line. 


43 


300 



1.0 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


Figure 6. Monthly averaged surface albedo (gray) and upwel ling shortwave radiative flux (dark) 
for SHEBA observed (solid) and corresponding M ERRA values (dashed) for October 1997 to 
September 1998, in W rrf 2 . “Tower” values are from downward-pointing pyranometer 
measurements, while “line albedo” values are from surface measurements along a 300m line. 
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Figure 8. Time series of monthly averaged near-surface station air temperature anomaly and 
corresponding MERRA valuesfor Barrow (left, 71°N, 157°W) and Jan Mayen (right, 71°N, 
9°W), in degrees C. 







Figure 9. Hourly (a) preci pi table water and (b.) liquid water path from SHEBA microwave 
radiometer and corresponding MERRA values, in mm. 
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Figure 10. Average annual time series for near surf ace station temperature and corresponding 
MERRA valuesfor (a) Amundsen-Scott (90°S), and (b.) Scott Base(78°S, 167°E), in degrees C. 
Bars indicate the standard deviation of monthly values over the period 1979-2005. 








